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ABSTRACT 



Field patterns are drawn for a plasma filled waveguide of circular 
cross section immersed in an axial magnetic field. The plasma is assumed 
to be cold and collisionless, and linearized field equations are used. 

y 

Besides the exact analysis the MHD approximation (low frequency) and 
the quasi-static approximation (slow wave) are used. Comparisons are 
made among field patterns determined by the various methods of analysis. 
Among the field patterns at different frequencies and wave numbers 
certain features are pointed out and comparisons are made. These field 
patterns include particle charge density and current density, total 
charge density and current density, the electric field and magnetic 
field. 
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FIELD PATTERNS FOR A PLASMA FILLED WAVEGUIDE 



I. INTRODUCTION 

The objective of this thesis has been to draw field patterns 
for a specific plasma filled waveguide in order that field patterns 
of plasma filled waveguides in general may better be visualized. At 
the beginning of Chapter I several assumptions are set forth. From 
these assumptions the necessary equations for the fields are derived. 
In Chapter II, the equations are organized in a manner used in arriv- 
ing at the field patterns and then specific regions of the ( iS-'U-J ) 
diagram are discussed. The ( 3 -u> ) diagram is divided into thq 

magnetohydrodynamic (MHD) region, the cyclotron resonance regiqn, 
the plasma resonance region, the region of cutoff fields (propaga- 
tion constant equals zero), and a region with the propagation 
constant intermediate between cutoff and resonance. 
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II. EQUATIONS GOVERNING FIELDS 

In this chapter the equations of the fields are derived. Eight prelim- 
inary assumptions are made concerning the plasma model and the typq of 
geometry. The first derivation is the dielectric tensor which comes from 
applying conservation of momentum principles to the particles. Under the 
heading of the determination of the fields the transverse fields are found 
in terms of the axial fields and the dispersion relation is derived. 

Under the heading for the determinantal equation boundary conditions are 
matched for fields having no variation in the angular direction. For this 
case the fields are written out in explicit form. Under the heading of 
cutoff fields a method of solution is given for finding the cutoff frequen- 
cies. A discussion of the MHD-regime is then made and finally the theory 
of the modified quasi-statics approximation is developed. 

2.1 Preliminary Assumptions 

The field plots of this thesis are for a very restricted situation; 
therefore the development of the equations will be restricted from the 
beginning and further restricted as the development of the equations 
progresses. A more general development which this thesis parallels 
closely is given in (l.) Waves in Anisotropic Plasmas , PP. (18-24), 

CH. 9, and CH. 10. The preliminary assumptions and restrictions are; 

(l) A linearized macroscopic description has been chosen for 
the analysis. This means that every time dependent field 
quantity that is determined is a small perturbation from 
an unperturbed state. 

(2. ) The system, waveguide and plasma is homogeneous in time 
and space. Fields of the form e^ 1x1 ^ will be assumed 
and the plasma characteristics do no change with position. 



6 . 



( 3 .) The geometry will be restricted to cylindrical coordinates 
with the radius of the waveguide a= . 05 m. 

(4.) The conducting surfaces of the waveguide are assumed to be 
perfect electric conductors, i.e. ? X E,! 0 at r : a. 

(5*) The plasma description is further simplified by assuming that 
no collisions occur and that the plasma is temperate. A 
temperate plasma is one which satisfies the inequalities that 
the particle induced velocity be much less than the random 
thermal velocity* which must be much less than the phase 
velocity of the wave. The lowest phase velocity which is 
considered is in the magnetohydrodynamic region where the 
phase velocity is as low as 1.75 X 10~ 2 times the velocity 
of light. For the electrons to have a thermal velocity 
less than one-tenth the phase velocity of the wave in the 
MHD region the temperature of the plasma must be less than 
10 1 * °K or in terms of electron volts the plasma must be 
less than a 1 ev. plasma. All the rest of the patterns 
have a phase velocity greater than one-tenth the phase 
velocity of light. 

*For a Maxwell velocity distribution the most probable 

velocity of a particle is given by uq, - ,1 2 Iff' where k 

1 m 

is Boltzman's constant, m is mass, and T is temperature. 

It is by definition the thermal velocity. 

(6.) The axial unperturbed magnetic field has a strength of .1023 
webers / m 2 . * 

(7.) The plasma is made up of electrons and ions, i.e., protons 

which have a mass of I .836 X 10 ^ times the mass of an electron. 
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( 8 .) The unperturbed density of each type of particles is 



*The values of the quantities in assumptions( 6 . ) and ( 8 .) were 
chosen as such for two reasons. The first is that a 
normalization of quantities which is used in the derivation 
of equations gives convenient numbers to work with. 

Secondly, the order of magnitude was chosen because they 
are thought to represent typical values for a plasma filled 
waveguide . 



2.2 Determination of Dielectric Tensor 

The equations governing the fields are Maxwell's equations and 
a conservation of momentum equation for each type of particles. 



where the k subscript stands for the type of particle. 

Since there is uniformity in the z direction, the fields will be 
broken into their transverse (T) and axial (z) components with^Bo - B 0 i z 
Solving equation ( 3 .) for gives 



In order to have the units stand out some normalization constants 
will be introduced and some quantities will be defined. 



1.808 X lOr 1 ? particles/m^* * 



\7 X - -J (jjjJq H, 

V X H = j o /£ 0 / J, 



Juan k vfc = e* (E,/ vk X Bo) 



( 1 .) 

( 2 .) 

(3.) 



J-W - ekBr> i z X 







£ 

a 



u 



ck “ 



e k ®o 



m k<% 



8 



( 6 .) 

(7.) 



V 



( 8 .) 



u 



u 




= u c 



+ Ur 



(90 



( 10 .) 



P ~ pe ' “ pi 

where c is the velocity of light in free space and a is the radius of 

the waveguide. u c ^ is the normalized cyclotron frequency, Up^ is the 

normalized plasma frequency, Up is the effective normalized plasma 

frequency, and u is the normalized frequency. Assumptions (6.), ( 7 .), 

and (8.) gives u ce = — 3 and Up » 4. 

Substituting the normalization constants into equations (4.) and 

( 5 .) and then writing the velocity in matrix form gives 

- ju 



v k =_!! 



m k w n 



I 



— u 



ck 



?_ ;“ 2 ck _ J _ L 



u ck 



I 



-j u 



u‘ 



.2 - u 2 



ck 



2 2 
u - u 



ck 



“j 



E (11.) 



U 

In order to determine J in terms of E some definitions are 

/\> 

introduced. 



A 


- N o ®k Vfc 


(12.) 


j 


= ie + 4i 


(13.) 


j 

/V 


= cr . e 

Aw /V 

A/ 


(14.) 


K 


= I + & 


(15.) 


/V 


~ j uo £ 



o 



jCT is determined from equations (ll. ), (12.), ( 13 . ) and (l4.). 
From equations (12.), ( 13 ,), and (l4.) 



* JL " N 0 ( e e^e + e A &) (l6.) 

Substituting equation (ll.) into equation (l6.) determines CT . The 

OJ thus found substituted into equation (15* ) determines K the dielectric 
tensor. Carrying out the above operations give 

9° 



K 



where 



K, I ~K X I 0 

- X - f - JL -I 



K x 

0 



r 



K.^ - 1 — 



u2 p (* 2 



K x I 0 

0 I K,. 



+ u ce u ci) 



( 17 .) 



(u 2 -u 2 ce ) fu 2 -u 2 ci ) 



= 



K = 



- j U U^ p (u ( 



ce 



u ci) 



( u2 *' u2 ce)(' u2 -' u2 ci) 
1 - /"P 



(?) 



(18.) 

(19.) 

( 20 .) 



Waves in Anisotropic Plasmas has the equations of this section on 
pages 18 - 23. Equation (ll.) is equation (2.8) on page 19, equation 
(17.) is equation (2.11) on page 20, and equations (l8. - 20.) are 
equations (2.21) on page 23. 

2.3 Determination of the Fields 

Substitution of equations (l4.) and (15») into equation (2.) 
leaves two vector equations with two vector unknowns. Equations (l.) 
and (2 . ) are rewritten along with two auxiliary equations found by taking 
the divergence of the two curl equations. These auxiliary equations are 
not independent of the first two, but nevertheless, are useful relations. 



V x 


E 


O £ 


(21.) 


V X 


H = 

A-' 


juu€ 0 K • E 


(22.) 


V ♦ 


K • E = 


0 


(23.) 


V 


0 H = 


0 


(24.) 



Again the procedure is to split the fields into transverse and 

axial components. Since uniformity exists in the z-direction, the axial 

- m 

dependence of the fields can be taken as e * Tbe operators become 
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v = v T - y t x (25.) 

4 • S. - K j.& + Kjc j*, x ET * K,, ^ H (26- ) 

After splitting the fields into their components, the trans- 
versely directed part of equations (21.) and (22.) can be solved 
for the transverse fields in terms of the axial fields. This is 
matrix ( 27 .). Using this matrix and the remaining parts of equations 
(l8 - 24) two coupled equations between E z and H z can be found in 
a straight forward manner. These are equations ( 36 .) and (37* )• The 
above procedure is carried out in more detail in Waves in Anisotropic Plasmas, 
pages 134 - 138 . Matrix ( 27 .) is the same as matrix ( 9 . 21) on page 135, 
and equations ( 36 .) and (37*) are the same as equations ( 9 . 30 ) and 
( 9 . 31 ) repsectively on page 136 . 



-Sr 




P 1 R j Q j S 




V I E z 






T | P 1 U 1 Q 




y T H z 




— 


1 1 ~ T 






r z x Ej 




1 

£> 

1 

CO 

S3 




r z x v T e 2 






-U 1 -Q 1 T Ip 

1 l 1 _ 




r x 7i\ 



where 



D = ( T + u K x ) + (u K x )‘ 

a 

P = a X 

P = -r(r 2 + u 2 ^ ) 

D 



( 28 .) 

(290 

(30.) 



(27.) 



R 

Q 

s 

T 

U 



J u3 ** 
D 



"HT 

\ £ o 



rAx 

D 

Ju (P 2 + u 2 ^) 

JP 2 u K x 
D 

-j u Kx (y 2 + U 2 
D 



\ € ° 




(31.) 

(32.) 

(33.) 

(34.) 

(35.) 
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( 36 .) 



The two coupled equations are 
V 2 E„ + c E„ - d H 



where 



+ f H z - g E, 



= (r^ + u^K.) _£* 

Jm 2 
a K, 



d = j u r K x 
a 2 K x 

^2 ^ ,.2 






N 



f = ( r 2 + u* K x + u 2 K v 2 



g - - j u P K|j 



a 2 K 



X 



fen 1 



>U 



(37.) 

(38.) 

(39.) 

(40.) 

(41.) 



Equations (36.) and (37.) can he combined to give two 
uncoupled fourth order equations. The equation for E® is 

2 



+ (c + f) Vt + (cf - dg)j E z = 0 



(42.) 



Plane wave solutions of the form exp (-jg*rr) satisfy equa- 
tions (36.), (37«) and (42.) where j? is a complex function of 
frequency only. The application of the identity (43.) to equation 
(42.) gives the dispersion relation for the medium equation (44.). 






e-J 




£ -x T 



(43.) 



p 1 * - (c + f) p 2 + (cf - dg) * 0 (44.) 

In order to satisfy the boundary conditions at the waveguide wall 
two two dimensional Helmholtz equations will be solved using 
equation (42.) to give the relation between the two P valves. 
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(V T 2 ■ -2 



Pi ) E zi " 0 



E z = E zj + e z2 

E z = h| E Z j + h2 E Z 2 



(45.) 

(46.) 

(47.) 



where is determined from equation ( 36 .) 

hj_ = c - p| 



(48.) 



The above equations can be found on pages (154.) and 
(155») in Waves in Anlstropic Plasmas . 

y 

2.4 DsterminaitaL Equation for Cylindrical Coordinates 

For cylindrical coordinates equation (45.) is 
^ 2 E 3i + 1 l E zi + 1 £ E zi + Pi 2 E zi = Q 

i 2 l IP 



(49.) 



b r 2 r ^ r 

Applying a separation of variables, equation (50.), to 

equation (49.) gives equations (51.) and ( 52 .). is the 
separation constant and is not to be confused with the m 
used for mass earlier. The m for mass will not be used in any 
of the following discussions. 



E zi - A i 


R(r) $(<N 




(50.) 


d 2 g 








d<£ 2 + 




0 


(51.) 


^£l. 2 + 

d (Pi^) 


p r dR(r) + 

d(pj.r) 


|(Pi r ) 2 -m^* 0 


(52.) 



In order to have the boundary conditions satisfied for all 
( 1 ) mj_ = m, and for the fields to be periodic with <p m must 
be an integer. Finally, R(r) - cfo(P ir) since there is no need 
for singular solutions in the waveguide. Therefore, the solution 
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sought is 



« 



E zi = A i J m ( p i r > eJm4) (53. ) 

For the rest of this section it is assumed that m = 0. 

For m ^ 0 see page 226. in Waves in Anisotropic Plasmas . The 
boundary condition E z = 0 at r = a gives 

A 2 J o (X 2^ = - A l J o< X l) (5M 

where 



Pi 



( 55 .) 



The boundary condition that E ^ = 0 at r = a gives the 

determinantal equation. E ^ is determined from matrix ( 27 .) and 
equations (46.), (47.) and ( 53* ) * The determinantal equation is 



where 



JiUi) 
bi J 0 Ui) 



bi = 1 



x 2 *^1 ( x 2 ) 

V J 0 < x a) 



( 56 .) 



k ± v 
a D 1 



( 57 .) 



Rewriting equation (44.) in its explicit form gives 



'•d21 




X 


- 


Q 


_ 1 
2 


(c + f) - 


(c - f) 2 +-4dg 


vt 


1 




2 









(58.) 

(590 



The simultaneous solution of equations ( 58 .), ( 58 .) and (59* )> 
for varying values of u determines the ( 6*u)) diagram. Fig. (l.) is 
the part of the diagram that is complete for m = 0. This diagram 
results from the work of Messrs. P. Serafim and P. Chorney. 

Using the following four definitions, the complex functions for 
the £ and Jtt fields for m = 0 are given by equations (64.) through 

( 69 .). 
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F = 
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\ J ° P. j) • 


•>“ <"nt) e - (rfl)} 


(60. ) 


F 2 " 


-A A J o »l) \ 


J Q £) . . -S 


Til! <«., 




Jo (^ 2 ) 


G 1 = 


a * r . 






J. 


Sr 




(62.) 


Go - 


a >*2 






2 


<5> r 




( 63 .) 


E = 
r 


2 

x <=* ' 


‘(P-U 2 K (1 )(p 2 + u 2 K_J+ U 2 Kj_ X ± 2 






i=l 


-P('r ,a " + u2 K j) 2 +(u 2 k^ 2 


(64.) 


E_. = 


2 

z - 

i=l 


'(PuK x ) 2 +('r 2 + u 2 Kj 2 K () - K L (p 2 +u 2 K x )x i 2 ‘ 




Hit* [(p2 + + (u a K I ) a J ’ 


E z - 


F + F 

1 2 




(66.) 



(65.) 





H 

z 




1 6 4 — 

i=l 



(PukJ 2 


+ (P 2 + 


2 , 2 

u Kj K, - K 


X (r 2 + U 2 K 


J ^ 


-juK x 


| fp 2 + 




1 : 


- 


“Ju [fax H 


+ 

CM 


u 2 K J +u 2 K x 2 - 


J^kXj 2 


(68. ) 



2 2 2 

(r + u kJ + (V k^.) 




C. 

I 

i=l 



(P 2 + 



U 2 K x ) K, 



K x X i 



ju pk^ 



(69. ) 



By inspection of equations (65.) and (660) give 




( 70 .) 
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A more general statement of equation ( 70. ) can be found from 



matrix (27.)° It applies to field patterns for which Vj E z ^ 
is parallel to Vx E Z g . By definition 

ti x T b - \ (71.) 

’ p i ( E J "5 (72-) 

Matrix (27«) then gives 

° Ejx = rb , i z X ^Q, + S Vj H ziJ 

(73.) 

v 

2 

^ ’ * “ T * ‘ X ( T E *i + P V T H zi ) 

>1 

H — 

Therefore a "* \ 

2.5 Cutoff Fields (Ps 0) 

By definition a cutoff wave is one whose propagation 
constant equals zero. When cutoff occurs the coefficients 
P, T* d, and g become zero and equations ( 36 .) and (37.) 

become uncoupled. 



go 



*b 



(75.) 



Therefore 

(Vj + p e 2 ) E z - 0 

2 2 0 

X e '(*.*) - u K || 

(Vj + P h 2 )H z . 0 




(76.) 

(77.) 

(78.) 



( 79 .) 



E - Waves : 

® s = A e ,J m (X e Cj 



h 2 . 0 

Em = 0 

* - ~ ja 
u 



x V T E z 



( 80 .) 

( 81 .) 

( 82 .) 

(83.) 



The boundary condition that E z =0 at r = a requires X g 
to be a root of the J m (X e ) Bessel function. The cutoff frequencies 
are given by 



H 



co 
Waves i 



H z - A h J m ( x h-5-) 



e ju(u> n t) 



E z = 0 



Hy= 0 



N 



X 2 h 






l k x 



v t h z + r 2 xv t h 2 



(84.) 

(85.) 

( 86 .) 

(87.) 

( 88 .) 



For the H=cutoffs the boundary conditions become more 
difficult; therefore m = 0 and m = 1 will be considered 

separately. For m = 0 V<j H z is in the T^, direction. The 

boundary condition E ( * ) =0atr = a requires 

T z x v t h z =0 ( 89 .) 
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Again the boundary condition requires to be a root of the 
J, (X, ) Bessel functions. The cutoff frequencies are given by 
equation (79*) which expands into equation (90.). Figures (2.) and 
( 3 .) are plots of equation ( 79 *) and nine of the solution points are 
indicated. 





6 

u co 


+ a^_ 


4 2 ^ 

u co + a 2 “co + a 3 ■ 0 


(90.) 


where 


a l 


2 

= u ce 


+ u 2 c i + 2u 2 p + X 2 


(91.) 




a 2 


CM 

11 


- u ce u ci) 2 + x h 2 ( u ce 2 + ^i +u p) 


(92.) 




a 3 


II 

&0 


u ce “ci ( u ce *ci " u2 p ) 


(93.) 



2.51 H - Cutoff (m = 1) 



V T H z = X n A h e ju(cx> n -t) 



L A h ' 






J i K 5) 

x hl 






] 



The boundary condition that E ( j,= 0 at r = a gives 



x h J 0 (x h ) - 1 - -J (95,) 

J (X h ) K i- 

The simultaneous solution of equation (79* ) which is 
repeated as equation ( 96 .) with equation (95 •) gives the cutoff 
frequencies for m * 1 





(96.) 



K X 

Figures (2.) and (4.) for low frequencies and figures ( 3 .) and 
( 5 -) for high frequencies demonstrate one method for solving for the 
cutoff frequencies. Equation (95») is plotted on figures (4.) and 
( 5 *)« The left hand side of equation ( 95 ») is plotted versus X^ . 



18. 



The right hand 6ide is plotted on the same axis to the same scale as the 
left hand side hut versus u . From these figures a function of 
(X h vs. u) is generated hy selecting an ordinant on the graph and 
reading the corresponding and u, then repeating with a new choice 
of ordinant. Equation ( 96 .) has been plotted on figures (2.) and (3»)» 

Also plotted on these two figures is the (Xh vs u) function generated 
from figures (4.) and (5«)« The dashed curves are for imaginary values of 
X^. Therefore, intersections of either solid lines or dashed lines 
simultaneous satisfy equations ( 95 * ) and ( 96 .) and thus are solutions 
from the H-cutoffs from m = 1. 



For finding field solutions in the cutoff region equations ( 83 .) 
and (88. ) are written in a more explicit form. Equation ( 83 . ) for 
m = 0 and m = 1 respectively is 



The explicit form of equation (88.) for m = 0 and m as 1 

respectively is 






( 100 .) 
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(4.) and ( 5 .). 

Waves in Anisotropic Plasmas considers the cutoff fields on 
pages 158 - 162 and page 227 * 

2.6 Low Frequency, MHD Regime 

By definition the MHD regime is the limit that occurs when the 
frequency is much less than both plasma frequency and the ion 
cyclotron frequency. In this limit 

K i — 1 - 

u ce u ci 

Kjj *«. -ju Up 2 (^e + u ci) 

— n 

^ce u ci 



(103. ) 
(104.) 



K, 



(c - f) - 

4 d g 



-f? 



1 + 






L X 

2 Up K* P 

a k! 



(105.) 



( 106 .) 



(107. ) 



In the limit of zero frequency ( constant), (K^' 



0), 



and (K 



II 



— Os5 ). A comparison of equations (106.) and (107*) shows 
» | 4 

c and Pg 



that | 4 dg| ; therefore from equations ( 58 .) 

and ( 59 . ) p^ ^ c and Pg ft; f. Only a slight coupling exists 
between equations ( 36 .) and (37* )• The MHO approximation considers the two 
equations to be uncoupled; 
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For the H-Wave, equation (37*) when (g = 0), the cutoff frequency 

is much greater than the ion cyclotron frequency, hence it is not in 
the MHD region. 

E - Wave : 

A/ 

V T 2 E z + p| E z * 0 (108.) 

x e 2 -(p c a ) 2 >(r > 2 + u 2 (109.) 

Equation (109.), the dispersion relation, is plotted in 
figure (6.). In this figure, the slope (Bd/u ) approximately 
equals 67, whereas a slope of 57*1 gives a phase velocity equal 
to the Alven velocity. The definitions of phase velocity and 
Alve'n velocity respectively are 
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V a 



^ce u ci' 



"p 1 
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(HI.) 



In the limit of zero frequency Q and T, equations (32.) and 
(34.), approach zero; and P and U, equations (30.) and (35* ), are 
given by 
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giving 
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(113.) 

(114.) 

(115.) 
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Equations (ll 4 .) and ( 115 .) were used in arriving at the field 
plots for the MHD regime rather than using actual values for the 
coefficients. The results seem questionable and will be discussed 
in the following chapter. 

2.7 Modified Quasi - Statics Approximation 

Near the resonances it is possible to have waves with arbitrarily 
slow phase velocities in the assumed plasma model. By definition the 
phase velocity of a wave is its frequency divided by its propagation 
constant. Quasi-statics is an approximation that takes advantage 
of the slow phase velocities. In this approximation the velocity of 
light is assumed to be infinite as compared to the phase velocity. 

From this it follows that the E^ field is approximately curl free. 
(This statement becomes more clear when one assumes plane wave 
solutions of the form exp.( — j Jk • £,) for Maxwell's curl equation 
for E. 

Av 

- J £ x £ a (116.) 



E 

r>j 






k 

X 

k c 

For finite H the right hand 






H 

/V 



( 117 .) 



side of equation ( 11? . ) approaches 
zero.) In Waves in Anisotropic Plasmas pages (162- 169) Professor 
Bers studies these resonances by taking limits of the exact solution. 

At plasma resonance he shows that the field approaches zero. 

For this case quasi-statics is a very good approximation. At cyclotron 
resonance he shows that (H z - 0 ) and JBy remains finite. For this 

case the dispersion relations differ in the two approximations. The 
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limit of the exact solution gives for the normalized p-number 



X 



ec 



JI M . + 



K 



2 u 2 K, 



(118.) 



while that found from quasi-statics is 



X 2 = p : 






(119.) 



Before proceeding with the derivation it should be noted 
that the field plots are for a modified quasi-statics. That is 
the field as determined by quasi-statics is used to determine 
an JR field, hence the name modified quasi-statics. 

Waves in Anisotropic Plasmas hns the derivation of the quasi- 
static field solutions on pages (240 - 248) and for the modified 
quasi-statics on pages (253-254). Also equation (ll8.) is the 
same as equation (9*219) page (l64.) in the book. 



2.7.1 Derivation of Field Equations 

For the quasi- static approximation equations (120.) and 
(121.) are assumed to describe the fields. 



V X E = 0 (120.) 

V. £ 0 K . E =0 (121.) 

Since the field is cure free a potential can be defined 

E = (122.) 
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Splitting the field into its components as previous proceedures 
followed give the following equations from equations (121.) and 
(122.) 



^ 7 0 + « 0 



(123.) 



P 2 = Y 2 JiL 



(124.) 




It ■' V t4 i 



(125.) 



B, 



z 



= Y <t> 



(126.) 



Figures (6.) and (7.) are a plot of equation (124.). The 
boundary condition that E z - 0 at r s a requires (X = pa) 
to be a root of the Bessel function. 

The H— field is determined from the E — field and Maxwell's 
equations. 



For equation (127.) a particular solution can be found which 
satisfies the right hand side and a homogeneous solution can be 
found which will be added to the particular solution in order that 
boundary conditions can be satisfied. 

Again the fields are split into their components and solutions 
of the form exp.(— £j.) are assumed. A few algebraic operations 
on equations (127.) and (128.) give 



^7 X H = J UJ £ K . E 



(127.) 



v. & - 0 



(128.) 
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y 2 $r 


(129.) 
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E z 




(130.) 


where Y s* -Ju 

T , 


*0 ' 
*0 


*&-** - 

V K " 


(131.) 


Y » 1 JL 

2 r 


I £o x 


(132.) 


Yo = -.1 u 

P 

The homogeneous part 


rzz k k (133.) 

fc ° Ax 

rsr 4 -k,, 

is found from equations ( 134 . ) and 


(135.) 

s' 








V X 4 - 


0 




(134.) 


V • A 


= 0 




(135.) 


Therefore Hh - — 


V(p 


(136.) 



v| y + r 2 v = o (137.) 

and the total H - field is 

H sH + H. (138.) 

Equation (138.) is complete when boundary conditions are 
specified. Since the boundary’s are perfect electric conductors, 
the arbitrary constant in the homogeneous part is found by 
lettering H r = 0 at r = a. 

After defining some symbols, the and JS fields are 
written in explicit form. 
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A( u> 6„ A t> = I v i G>H*> - o,«-S-)].3-o 



(139.) 



Z = X 



TkT 



K, 



C h = 



Y 2 J 0 (x) 



lo(6„) - 



iiOV 

6 , 



m ■ O: 



E = A(u,6 n# t) 



Kz A(u.fl^t) 



j^X 



• J l<* i> T r + V* -I) 

~ IJTT) r . 



JV 

4, 



+ fxr 5 j.d - r ) - 4444 i-o(e„- r ) 1 r. 



m = 1: 



E = A(u,0.<p,t) 
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a 

x J i (X r ) -.1 

-f ^rt- r * - ; J *. 



° v 



(140.) 



(141.) 



(142.) 
(143 . ) 



(144.) 



£ = A(u,e^,t) r ± r + 4 * h z rj 



where h r = -j 
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III. FIELD PLOTS 



The section describing the conventions has the equations which 
describe the particle motions listed and a discussion of the two 
conventions used in generating the field patterns. The two conventions 
are (l.) polarization ellipses at various locations to describe the 
field strength and a vector to describe its direction and (2.) spacing 
of lines to determine the strength of the field and direction of the 
lines to describe the direction of the field. All the rest of the 
sections consider the plots for various regions of the (B-uj) diagram. 
Comparisons are made between various plots and the interesting 
characteristics of certain field patterns are pointed out. 

3.1 Equations and Conventions Used 

The objective of this thesis has been to use those equations 
which were developed in Chapter I in order to arrive at patterns 
of the fields for selected points of the (B-VD) diagram. Figure (l.) 
has indicated on it those points for which field patterns have been 
drawn. In order to get the equations for the field pattens and 
II were determined from those equations developed in chapter one in the 
various sections. The rest of the field patterns were calculated by 
using the following equations. From equations (l4.) and ( 15 .) 
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Equations (8.), (ll.), and (12.) give 
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0 ! 0 i 
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J = J. + J 



(151.) 



t \ = “77 



UUJ 



? = u 



UJ 



n 



n 






V • J 

r-' 



(152.) 

(153.) 



Taking the curl of equation (2.) and equation (153.) give 



■P = V . € 0 E 



(154.) 



The constants for the plasma are given by 



u pe _ ^ 



= 



16 



1836 



- ^ = 4(155.) 

.1 * 



u ce " ~ 3 u ci = -I,. 

1836.1 



(156.) 



Equations (151.) and (153.) were used to check the correctness 
of the answers. Equation (151.) was vised to check that sum of the 
current densities as determined from equation (150.) actually 
equalled the current density as found from equation (l49<). 

Equation (153.) was a second method for calculating the total 
charge density and was used to check the correctness of the answers 
as found by equation (154). This check is equivalent to equation 
(23.) 7 . K . JE^ = 0. Finally, since the computations were 

performed by slide rule, no more than three significant figures 
can be expected. 

No patterns were drawn for the separate particles except for 
frequencies below (u - l). This is because nearly all the current 

is from movement of electrons at high frequencies. Therefore, the 
pattern for the electrons is the same as the total current. The pattern 
for the ions is nearly the same as that for the E - field. 
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2 2 

The justification of this statement comes from assuming u u 



and looking at equation ( 150 .)* 




\ u / a 




The most noticable case for disagreement between the patterns for 
E^and Jj. would be for ^Erp to be linearly polarized. This would 
produce elliptical polarization for^j^, with the ratio of major 



the field pattern and the current density pattern for the ions is 

approximately the same as that for the E field. 

r- 

For the field patterns two sets of conventions were used to 
give two sets of field patterns. The first convention has ellipses 
of polarization at various locations to show the relative strength 
of the fields. In each ellipse a line is drawn from the center to 
the outside which indicates the strength and direction of 'the field 
at time t = o when the constant of amplitude for the fields is 
taken as a positive real number. As time increases the arrow 
indicates the direction in which the vector rotates. 

The three components of any field are drawn to the same scale 
unless it is noted otherwise by (X Number) on the pattern. The 
particle current and the total current are drawn to the same scale. 
Finally, the scale of any particular or ^ ) was chosen to 

give convenient sized ellipses of polarization. 



to minor axis equal to 1.6 X 10”^. Therefore, it is concluded 
that the term ( u c j_ X E) has negligible affect on 
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Only two of these views are necessary since the third can be 
drawn from the first two. The (r - <()) view is taken looking down 

the waveguide axis ( / z -direction) at z = 0; the ( r - z) view 

is taken looking across the waveguide into the minus ^-direction 
at <f) = 0. The ( ( <3>-z) view actually is generated from the preceeding 
views. It corresponds to looking into the minus r-direction at 
4> s 0 and moving the polarization ellipse up or down a given dis- 
tance r. depending whether the ellipse is at plus r or minus r 

respectively when 4> - 0. 

The second convention indicates the strength of the field 
by the spacing of lines. This convention has been applied to all 
of the transverse fields for m “ 1 and to the radial fields 

for m * 0. The ^-directed component is divergence free for 
m - 0 and its fields are readily visualized from the patterns 

of the first convention. The z -dependence of the patterns was not 
accomplished due to lack of time but can be sketched easily from the 
patterns of the first convention. For an increase in z the vector 
rotates in the direction of the arrow. When z - TT the 
vector has rotated by 90° in space. 

The second convention illustrates two features. It clearly 
illustrates the divergence of the fields, especially at cutoff. 

Laying the pattern for the jS- field or for the Sj field rotated 

by 90° in space over the charge density contours shows that where 

the density is the greatest the number of lines starting or ending 

per unit area is also the greatest. The second feature is that the total 

transverse pattern rotates in the counterclockwise direction for m ~ 1. 
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Before proceeding with the discussions of the individual regions on 
the (0-w) diagram, it should be noted that the assumed model of the plasma 
may break down at the boundary. Inherent in the model is the assumption 
of specular reflection. Specular reflection accounts for the fact that 
particles cannot flow through the walls of the waveguide as the field 
patterns indicate. 

3.2 Low Frequency, MHD Regime 

For this region, four sets of fields were calculated, two by 
the MHD approximation and two by the quasi-static approximation. 

Such good agreement occurred between the two approximations that only 
two sets of patterns were drawn, one for m - 0 and one for m - 1. 

The largest variation between the two approximations was a factor of two 
for and ^° 

For this region and exact analysis was not possible because of 
failure to find a solution point on the (8-ul) diagram. A comparison 
of the coefficients which were used with those from an exact analysis 
using the same frequency and propagation constant will indicate the 
degree of the error. For an assumed u = 17 X 10 _ 5 and/3 n » 1.2 X 

10”^ one of the p numbers^ equation (58.), is very close to c, equation 
(38.), but for the exact analysis it is not necessarily a root of the 
Bessel function since there are two p numbers. Coefficients which would 
have been used for the approximate analysis to find the transverse fields 
are 



U = 1.53 



T 2 0.0 



P = j 12.3 



Q - 0.0 



(158.) 

(1590 

(l60.) 

(161.) 
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while the coefficients from an exact analysis are 



p = 


j 12.9 


(162.) 


Q = 


-2.6 


(163.) 


T 


0.5 


( 164 .) 


U = 


Jl.63 


( 165 .) 



From a comparison of the coefficients and knowing that there 
are two p-numhers each of which are different from the approximations 
raises an element of doubt as to the value of the approximations. 

Assuming the approximations are a valid representation of the 
fields, it can then be concluded that for both m ■ 0 and m = 1 
the fields lie in the transverse plane except for jJ and Jg . The 
strong axial electron current is a result of (Kyj — ►- °«) as (ui— *-0). 
Another interesting feature of this pattern is that the ions move in the 
transverse plane in such a fashion that the divergence of the total 
current is negligible compared to the divergence of either particle . 
current, i.e., the particle charge density is approximately 3 X 103 times 
greater than the total charge density. 

3.3 Electron Cyclotron Resonance 

Three sets of field patterns, two for m 2 1 and one for 

m - 0, were drawn for this region from the quasi-statics approxima- 

tion. Here larger p-numbers than those used for other regions were 
assumed for two reasons. The first one was to develop some patterns 
with more variations in the radial direction. The second and most 
important reason, is that the dispersion relation as determined by quasi- 
statics approaches the true dispersion relation only for large p numbers. 
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This can be seen by a comparison of equations (ll8.) and (119.). The 
model for the quasi-statics approximation assumed the velocity of light 
to be infinite compared to the phase velocity of the wave ? in this case 
the phase velocity in the transverse direction, hence large p numbers. 

Two points seem to be indicated from these patterns. The first is 
that some of the lines for Jp and Hr enclose upon themselves. The 
for other regions is completely encircling, cutoff and MHD, or no 
encircling lines as plasma resonance. The second point is that a 
comparison of the fields at electron cyclotron resonance with those 
at plasma resonance shows that the H-field for electron cyclotron 

/N/ 

resonance is about an order of magnitude greater than the H, - field 
at plasma resonance when the two fields have similar propagation 
constants and similar values for . 

3-4 Plasma Resonance 

This is an interesting region because several comparisons can 
be made. As a first choice the exact analysis and the quasi- 
static analysis at u s 4.08 will be made. Here it can be seen 
that there is very good agreement between the two patterns and field 
solutions. The difference which occurs is that the E - field becomes 
slightly elliptically polarized in the transverse plane for the exact 
analysis and the strength of the z -directed field as found by the 
exact method is about twice that from the quasi-static method. 

The second interesting comparison is that at similar propagation 

constants the plasma resonance canes closer to being an ^E-wave than 
1 < 

cyclotron resonance as was discussed in the preceeding section. Also 

i* 

E z for plasma resonance is about twice E z at cyclotron resonance when 
the two have similar magnitudes for 
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The third interesting comparison is to compare the patterns of those 
fields which have solution points that lie on the curve on the (6-ud) 
diagram that originates at u Z 3 . 94 . See Figure (l. ) . As the propagation 
constant increases the z component of the E_ and J fields grow while the 
H_ - field decreases. The transverse JH - field grows from nothing 
at cutoff to strong radial then to almost circular while the E - field 

/V 

starts frcm almost circular and decreases to almost radial. The ^ field 
is very similar to the Rg field except that it has a strong radial field 
at cutoff. 

The fourth interesting point is that J and H have a polarization 
in the transverse plane which is opposite to that of the field. 

The E vector rotates in the counterclockwise direction while J and H 
rotate in the clockwise direction. 

3.5 Cutoff Fields (pa 0 ) 

. 3 . 5.1 E - Waves 

All that seems apparent for these cutoffs is that the R field is 
the only one with transverse components and the charge density is zero. 

They are similar to the cutoffs for E - waves in an empty waveguide. 

They differ in that the cutoff .frequency is higher than that for the 
empty waveguide and that the empty waveguide has no particles, hence 
no current density. For the empty waveguide (u<:o- Xe ) of equation 
(8.0) and is a root of the Bessel function. 
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3.5*2 H - Waves 



m = Os 

Three field patterns were drawn for m = 0 . For the cutoff 

at u = .0430 the ion and electron charge density nearly cancel 

each other. From the patterns of and the particle paths 

can be drawn for low excitation, i.e., localized motion of particles. 




(a.) The small ellipse is the path of the ion. 

(b.) The large ellipse is the path of the electron. 

For u = 3.94 the interesting feature is that a nearly 

circularly polarized - field produces a strong radially polarized 
- field. The reason for this is that u ce E r approximately equals u E<^ 
in equation (150. ) hence (j^ ^0). 

The cutoff frequency u = 6.34 has nothing that particularly 

stands out. A comparison of the three patterns shows that as frequency 
is increased, the E-field goes from nearly radial, to almost circular 

with the stronger field being r-directed and then to almost circular 

with the strongest field being 0 directed. For m = 0 the H wave 

at cutoff in an empty waveguide has the field entirely in the 
<t> direction; therefore none of the patterns are like those of the 
empty waveguide. The higher the frequency the more nearly the modes 
look like empty waveguide modes. 
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m = 1: 



For m * 1 there are four sets of field patterns. A comparison 

of the patterns for U = .04326 and m = 1 with U = .0430 and m ~ 0 

shows that for both sets the charge density of the electrons tends to cancel 
that of the ions, that is, p e = — 3*75^> and p^ ■ 4.75p. The 

particle paths for m = 1 look like 




(a.) The smaller ellipse is the path of the ion. 

(b.) The large ellipse is the path of the electron. 

In order to visualize the way the electrons tend to cancel 
the charge density of the ions the second convention is the better. 

A comparison of u r 3*37 for m * 1 with u s 3.94 

for m - 0 shows that both have nearly circular polarization 

ellipses across the cross section of the waveguide, and both have 
the same direction of polarization. 

The last two sets of patterns are for imaginary p numbers. 

The fields for U • 5»07 show how large imaginary p numbers produce 

very strong fields near r ~ a. This is because a modified Bessel 
function grows rapidly as its argument grows. The fields for 
u - 5*77 are for p numbers that are practically equal to zero. 
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3.6 Patterns Between Cutoff and Resonance 



Two sets of patterns have been categorized under this heading, 
one for u ~ 4.773 and one for u a 1.00. For the patterns at 

u Z 4.773 one of the normalized p- numbers is 8.65 and the other 
is 2.40. The large p -number causes three variations in the field 
pattern as r goes from zero to a. For this pattern it is neces- 
sary to use the patterns from the second con'vention because the large 
number of variations renders the first set of patterns practically 
useless. 

The special feature the patterns at u - 1 have is that this 
is the only set which was drawn for complex p — numbers . For it 
a new complex constant of amplitude was chosen in order that the 
fields would be in either the (r or <t> )-directions at time 
t s 0. Like the cutoff fields with imaginary p - numbers the 
fields are strongest. near the outside of the waveguide. 
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MHD: 

E-Wave (1. ) pp. 169-172 
m = 0 

u = 17 x 10“^ 

3 n = l.llp. x 10 ~ 2 
p = 2.i|.05/ a 

A(*,t) = v J(17xl0- S ui n t-l.llOxl0- 2 1) 

Complex Fields? 

E * A(z,t) [-J795 J x (pr)t r + J Q (pr)? 2 ] 

H = A(z,t) [-J103.8 J^prjt^] 

J = A(z,t) [23.7 (pr)t r + J2.38 J^pr)*^ 
-jif.99xlO +3 J 0 (pr)t 2 ] 

J e = A(z,t) [+1.273xlO “ 2 J 1 (pr)t r -j2.25xl0 2 (pr 
”04»99xlO +3 J o (pr)f 2 ] 

= A(z,t) [23.7 (pr )? p + J2.27xl0 2 ^(pr)?^ 

-j 2. 72x10"^ J Q (pr)t 2 ] 

P = A(z,t) [-J3.38xlO“ 7 J 0 (pr)] 

P 0 = A(z,t) [-jl.ll8xlO“ 3 J Q (pr)] 

Pi = A(z,t) [ +jl. Il8xl0’ 3 J 0 (pr)] 
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%iasi-Static: (1. ) pp. 241-246, 253-254 

TO — 0 

u = 17x10“^ 

/ 3 n = 5 . 87xlo“ 3 , / 2 -= .1174 

p = 2.405/a 

A(z,t) = Ai e^ 1 7xlO" 5 UJ n t-5.87xlO’ 3 2 /a) 

Complex Fields: 

I = A(x,t) ( - J4l0 J x (pr )^ r + J Q (pr)t z ] 

H = A(z,t) [ (-6. 44x10 ^(prJ-3. 35x10”-^ J)? r 

-J103.9 ^(pr)?^ +j(2. 64xlO” 2 J Q (pr) 

-I.l 42 xl 0 ” 2 )l ] 
z 

J = A(z,t) [12.20 J^pr )T r + jl .230 (pr)lt^ 
-j4*99xl0 3 J Q (pr)t z ] 

1 Q - A(z,t) [6.56xl0” 3 (pr - jl. I59xl0 2 J 1 (pr)Tj ) 
-j4* 99x10 +3 J o (pr)t z ] 

J ± = A(z,t) [12.20 J x (pr )"? r + jl,171xl0 2 J-^pr)?^ 

-j2. 272x10"^ J Q (pr)t z ] 
p = A(z,t) [-jl. 746 x10“7 J q (pr ) ] 
p e = A(z,t) [-J5. 75x10’^ J 0 (pr)] 

Pi * A(z,t) [+J5. 75x10"^- J 0 (pr)] 



Note: Field Patterns Drawn For MHD Only 
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f ' (J>*0,Z=0 

L - > -t 
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.1 




u>i= -W 

P=P« 

P e = 3. 3 |X|^p n 

ft = -3.3W lo 3 P 
h 
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MHD: 

E-Wave (1. ) pp. 169-172 
m = 1 

u = 17 x 10 “^ 

= 1.353x10” 2 

p = 3 . 832 /a 

A(4»z,t) = v Ja7xlO-5w n t-1.353xlO-Va) eJ+ 



Complex Fields: 

2 Jt (pr ) . 

E = A(<>,z,t) [j5.93xlO^(J 0 (pr) - )t r 

2 J l <P r ) ^ 

- 5.93X10 2 -3~r- + ^(pr)^] 

J, (pr) J T (pr) 

H — A(^>,z,t) [65.1 + j65. l(J 0 (pr ) - pp ^3 

J-. (P** ) 

J = A(4>,z,t) [-17.63 J Q (pr) + 15.85 ■■ p — n r 

J l (P r 3 .* 

-3(1.779 J Q (pr) + 15.85 - pp— 

-j4»99xl0 3 J^pr)"^] 

-7 p J i (P r ) 

J e = A(^z,t) [-(9.49x10 3 J Q (pr) - 1.672x10^ ■ - 1 — -)t p 

2 J-iCpr) . 

+jl.672xl0 2 (J o (pr) - 




- j4»99xl0 3 J-^pr)^) 



A(<|>,z,t) [(-17.62 J 0 (pr) - 1.519xl0 2 

-j (1« 697 x 10 2 J 0 (pr) - 1. 5l9xl0 2 
-j2. 72 x 10 *’^' J-^pr)^] 



J x (pr* ) 

pr 

J x (pr ) 
pr 





p = A(4>,z,t) [«j4.02xlO~ 7 J-^pr)] 
p e = A(<l>, 2 ,t) [ -jl. 32ipcl0 “ 3 J-^pr)] 
Pl = A({>,z,t) [+jl.325xl0 “ 3 J x (pr )] 
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Quasi-Static: (1.) pp. 241-246, 253-254 

m = 1 

u = 17x10“^ 

B n = 9. 36x10 " 3 £= .1872 
P = 3.83/a 

A(*, 2 ,t) = ^.JUTxio^Wnt-go&cio -3 */*)^ 

Corqplex Fields: 

J-, (pr ) . J n (pr ) 

E = A(<|>,z,t) [jif09(J o (pr) — — )t r - lf09 — 

+ J^pr)^] 

- , MP 1 *) -* J i <P r ) ^ 

H = A(M,t) [65.1 “~F- % * J65.1(J 0 (pr) - -±— )t^ 

+jl.652xl0" 2 J, (pr)t ] 

Z J-. (pr) 

1 = A(*,z,t) [(-12.19 J 0 (pr) + 10.95 )t r 

-Mpr) 

-J (1.228 J Q (pr) +10.95 

-j4*99xl0 3 J x (pr )t 2 ] 

o ? J n (pr ) 

J a = A(^>,z,t) [-(6.5lpcl0*’ 3 J (pr) - 1.156x10 — — )t 

2 J (pr) ■* 

+jl. 156x10* (J 0 (pr) - 

-j4-99xl0 3 J x (pr y? z ) 

p J. (pr ) 

J A = A(<j>,z,t) [(-12.19 J Q (pr) - 1.048x10* -±—)l p 

P P Jn (pr ) 

+j (-1.169x10* J Q (pr) + 1.048x10* 

-2. 72x10 J x (pr )? z ] 

p = A( <|>, z, t ) [-j2. 77xlO~ 7 ^(prJJ 
p Q = A(<J>,z,t) [-j9. 16xl0“^ J x (pr)] 
p ± = A( {>, z, t ) [j9. 16x10 J x (pr )] 



Note: Field Patterns Drawn For MHD Only 
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electron CYCLOTRON RESONANCE: 

Quasi -Static Analysis 
m = 0 
u = 2,ij. 

f 3 = 10.08 ; < 3 « 201.6 

n 

p =5|S 

A(z,t) = V J( 2 -4"nt-10.08 f) 

Complex Fields: 

E = A(z,t) M.54 7 ^(prjTp + J 0 (pr)t z ] 

H = 10~ 3 xA(z,t) f (1. 6lp. (pr ) + 1.937x10“^ ( ,3r))? r 
. -J2.05 Ji(pr)t^ - j (.898 J 0 (pr) 

+ 1.937*10 I o ( 6r))t z ] 

J = A(z,t) [.3i|4 J x (pr )? r - j.lj-30 J 1 (pr )3 l^ 

-J.354 J 0 (Pi*)t 2 ] 

p =A(z,t) [-j2.23xl0 “ 9 J 0 (pr)] 
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ELECTRON CYCLOTRON RESONANCE: 

Quasi -Static Analysis 
n = 1 

u = 2.4 

Q n = 7.01 ; / 3 = 1 ^ 0.2 

p = 2^1 

A(<{>,z,t) = A 1 e^ ^■ u; n t “7. °i 5)^ 

Complex Fields: 

J, (pr ) J- (pr ) 

E = A(4>,z,t) [j.547(J 0 (pr) - )t r - -547 ^ 

+ J^prjtj 

.o , *Mpr) 

H = 10 j A(4>,z,t) [(-2.36 J Q (pr) + 5.31 - 1 — — 

- I, ( 0r ) 

- 6 . 41 pao“ 3 (i o ( 0 r) - ~ -g- r - - )? r 

J-,(pr) . , M er U 

+ j (2.95 J Q (pr) - 5.31 -±£5 6.44x10 3 )tj 

-j (1.292 J^pr) - 6.44 x 1°“ 3 I 1 (/3r))^ z ] 

J-, (P** ) 

J = A( $>, z, t ) [(-.34^- J 0 (pr) + .744 -—T— 

J-i (pr) 

+M.430 J Q (pr) - .744 
-j.1937 ^(pr)^] 

p = A( <j>, z, t ) [-J1.611X10" 9 (pr )] 
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electron CYCLOTRON RESONANCE: 
Quasi-Static Analysis 



u =2.18 

/3 n = 9.93 &= 198.6 

A(M.t> = V i< 2 - l8uJ nt-9.93 S)^ 

Conplex Fields: 

B = A(4>,z,t) [j.707(J o (pr) - -=^5 — )t r - .707 ^ 

+ 

- J- (pr ) 

H = 10”^xA( <J>, z, t ) t (-1.425 J Q (pr) + 3.38 ■■■ +— - 

). I n (6r) 

+ 1.794x10' % 0 (i3r) - -i -g- — )t p 

+j (1.962 J 0 (pr) - 3.38 ■ - - --- + 1.794x10 ^ 

-j (1.003 J x (pr ) + 1.794x10'^ I 1 ( (Br ) )t z } 

_ J-,(pr) 

J = A(<|>,z,t) [(-.308 J 0 (pr) + .732 ■ y-" )^ 

, , J l(P r ) 

+j(.424 J 0 (pr) - .732 

-j.390 J-^pr)^] 

P = A(<j>,z,t) [-j2.63xl0* 9 J-^pr)] 
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PLASMA RESONANCE? (1. ) pp. 162-I6if 
Exact Analysis 
m = 0 



u - if. 08 

G n = 10 ' 53 




p 2 * 2 

A( z, t ) - ^^4.08^.10.53 I) 



Conplex Fields: 

E * A( z, t ) [-J2.09xl0“ 1 (J 1 (p 1 r) - 6. 35*10‘ 5 1 ) )t p 
+ 3.9ipclO” 2 (J 1 (p 1 r) - 2. 4.9x10“^- I^^r))^ 

+ (J 0 ( P]L r) + l.l|27xl0 “ 5 I 0 (^ 2 r ))^ z 3 

H = A(z,t) [-2. 71x10”^ (J^p -j^r) - 2.1|.7xlO' , ‘^ 
+jl.722xlO“ i ^(J 1 (p 1 r) + 3.25xlO " 6 
+j6.25xl0“ 5 (J o ( Pl r) - 9.89x10"^ W^z 1 

J = A(2,t) [-8.l5^10“ 2 (J 1 (p 1 r) - 3.36x10“^ I 1 (Q 2 r))t p 
+j5.17xl0’ 2 (J 1 (p 1 p) + 3.9 x10" 7 
=.j2.08xl0“ 1 (J o (p 1 r) + 1.427xlO " 5 I o (Q 2 r))t z ] 

p = A( z » t ) [-jxlO" 9 (1.9^8 ^(p^) + 3. 76xlO” 6 I 0 (Q 2 r ))] 
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PLASMA RESONANCES 

Quasi-Static Analysis 
m = 0 
u =1^.08 

/3 = 164.6 

j(l|..08u> n t-8.23 f) 



= 8 a 23 r 



A(z,t) = A^e' 



Complex Fields : 

E = A(z,t) [ -j 2. 92x10 ’ 1 J x (pr)t r + J Q (pr)? r ] 

H = A(z,t) [- 5 . 69xlO’^‘(J 1 (pr) - 1. 045xl0” 3 1^ ( ) )”t r 

+j4* 20 x 10“^ J^(pr)^ 

+j6.92xlO“^(J (pr) - 8.6lxl0" 3 I (6 r))?l 

Q O Z 

J = A(x,t) [- 1. 322x10 -1 J^prjt + j9.72xl0~ 2 J-^prJlty 

-3. 208 J 0 (pr)t z ] 

p = A(z,t) [-jl.£ 8 xlO ~ 9 J Q (pr ) ] 
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PLASMA RESONANCE - PAST WAVE: 



Exact Analysis 



ra = 0 



u = 4*154 



£ n = 4.00 

Pl = 3.l49/a 



P 2 = jl. 874/a = j q 2 

, (4. 154 tu t-4 

A(z,t) = n 




Complex Fields: 

E = A(z,t) [j(-.686 J-^pp) + .0927 I 1 (q 2 r ))t r 

+ (.393 J^P^r) - .0782 I^q^r))*^ 

+ ( J 0 < p l r) + *14^4 I 0 (92 r ))^ z ] 

H = 10" 3 xA( z, t ) [ -1. 006 J x (p^p ) + . 199 7 I x ( q 2 r ) )* p 
+ j ( * 252 J^r) + .0620 I-^p))^ 
+j(.795 J 0 ( p l r > “ • °937 I 0 (q 2 r))t z ] 

J = A(z,t) (-.1717 ^(P^) + .01545 ^(qgr))^ 

+K.0437 j^p) + .00475 i 1 (q 2 I, ))‘ J < j ) 
-j(.2o4 j 0 ( P;l p) + .0300 i 0 (q 2 r))* z ] 

p = A(z,t) [-jlO" 9 (1.090 J q ( Pl r ) + .0728 I 0 (q 2 r))] 
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PLASMA RESONANCE: 



Quasi -Static Analysis 



m = 1 
u = if. 26 

B n =9.62 ; 8 = 192.4 



D = 
p a 

A(4>,z,t) 



. (4» 26 uj t-9. 62 

v 



Complex Fields: 

_ J-,(pr) J, (pr) 

E = A($,z,t) [j.398(J 0 (pr) - ■■ p y -')? r - .398 1 pr 

♦ J^prlTj,] 

? J.. (pr ) 

H = 10^xA(4>,z,t) ((.497 J 0 (pr) - .847 ■ - p ' r - - ■ 

-4 „ M Qr > + 

+ 1.129x10 ^I 0 (8r) - ■■-■ fly " ])*, 

. Ji (P r ) . I,(8r) 

-j (.350 J 0 (pr) - .847 1.129x10 ^ ■■ ■ 



+j (.1977 ^(pr)- 1.129x10^ I 1 (6r))t z ] 

_ J-,(pr) . 

J = A(*,z,t) ((.1571 J 0 (pr) - .268 -- 1 — -)t r 

J-,(pr) 

-3 (.1108 J Q (pr) - .268 )t^ 

-J.1992 J!(pr)t z ] 

p = A(4>,z,t) (-J1.970X10 " 9 J x (pr ) ] 
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CUTOPP FIELDS: 

E-Wave 
m = 0 

u = i|_. 67 

P = 2^ 

F a 

A(t) = A 1 e^- 67a ’ nt 

Complex Fields: 

E = A( t) [ J 0 (pr)^ z 3 

H = A( t) [j.1367 J^pr)^] 

J = A( t) [-J.1818 J 0 (pr)^ z ] 

J Q « A(t) [-j. 1818 J 0 (pr)* z 3 
1 ± = A(t) [-j9. 89x10’^ J 0 (pr)t z ] 
P - 0 
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CUTOFF FIELDS: 

E-Wave 
m = 1 

U = 5.5k 

p =M2 

y a 

A(*,t) = .J* 

Coup lex Fields: 

E * A( 4», t ) [J, (pr)l j 

5 = 1.833x10 3 A(4>,t) I- % - j(J 0 (pr) - -1—)^} 

J = A($,t) C-J.1534 Ji(pr)t z ] 

P » 0 



E 4 J 




tl r _(K>,z o 
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CUTOFF FIELDS: 

H-Wave 
m = 0 

p = Ml 

y a 

Conplex Fields: 

(1.) u col = .01^30 A(t) = 

E = A(t) J x (pr) [ 2 . 72 x 10 2 t p - J4-23 3M 
H = A(t) J 0 (pr)t z 
J = A(t) J x (pr) l-j. 619 t r + 77.0 
7 q = A(t) J-^ (pr ) [+J2.30 t r - 77.0 t^] 

= A(t) (pr) [-J2.92 t p + 6. 56x10 “ 2 t^] 
P = A(t) J 0 (pr) [1, 83xlO" 7 ] 



(2.) u co2 = 3 * 9 ^ A(t) = 

E = A( t) J x (pr ) [-4.98 t r - J3.88 t 

H = A( t ) J 0 (pr)t 2 
J = A( t) J-^pr) [jl04.0 t p - 4.48 
P = A( t) J Q (pr) [-3.37xlO~ 7 ] 

(3.) u cq3 = 6.34 A(t) = 

E = A( t) J x (pr ) [310 t r - j624 
5 = A( t) J Q (pr )t 2 

J = A(t) J x (pr) [-J104.6 t v - 133.2 
p = A(t) J Q (pr) [2. 10xl0“ 7 ] 



Lfe.0430 E 





r_ r . «--0,Z-0 C 
L “*V» l <f> 



f 



/ 100 



1 



r 



•r 





102 



U=043 0 E r 4J r 4J ffr 4j ir 




103 




104 



U=.04 3 0 

Co/ U 




tr *=«'0r 



X <?5 

h 






i 

i 




t •=£? 

p --e« 

^-- -&7ir h 

?■---> *-. i\ 



x^s 



105 
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CUTOFF FIELDS: 
H-Wave 
m = 1 

u = „0ij.326 

— A-. © 



j.0433W n t 



Complex Fields: 

Jl(pr) . 

E = A([>,t) [-(272 J 0 (pP) " 27S p r )t r 
+j(l|..17 J 0 (pr) - 275 pr 

H = A( ^ ,t ) [ 

J = A( <J> ,t ) [ j (.615 J 0 (P1*) + 76.5 p" ' )T r 

J-, (pr 5 ^ , 

“(76.9 J Q (pr) - 76.5 p y- )ty 

. MP**) ^ 

J e = A(*,t) [ - j (2. 29 J Q (pr) - 79.4 pr —'^r 
-(76.9 J Q (pr) - 79.4-^““)^] 

, J i ( pr > .* 

J ± = A(4>,t) £ j (2.90 J Q (pr) - 2.84 ■ ' ~ - )T p 

J,(pr) . 

-(6.52x10“^ J Q (pr) + 2.84 pr - -)^] 




t!: r $o.z=o £ 



T 



it 



i* 

H 



x 




II 

I* 







2 



i 

1 

\ 

i 

H 



f 

! 

f 



k 



l 



N> 



115 



J T 
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CUTOFF FIELDS: 
H-Wave 
m = 1 



u = 3.37 

- 2.42 
n — ■■ * 77t 

H a 

A($,t) = A^ 3 ’ 37 "*^ 



Complex Fields: 

J l(pr) . 

E = A(<t»,t) [548 J 0 (P r ) “ 2 3* 1 “pF~ )T r 



J-,(pi*) ^ 

+j(52lf J Q (pr) + 23.1 p r — - )^3 

H = A( <J>,t ) J 1 (pr)t z 

J 1 (pr 5 x t 

J = A( <{> , t ) M(98.0 J Q (pr) - 52.9 pr ™ A r 

J 1 (pr > x* 

+(144.9 J 0 (pr) + 52.9 pr *"" * * 4 

p = A(4,t) [-2. 35x10 ~ 7 J x (pr)] 
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CUTOFF FIELDS? 

H«*Wave 
ra = 1 
u = 5.07 

A( <(> > t ) = A 1 e j5 *° 7UJnt e M 
Complex Fields: 

E = A($,t) £(1.985 x 10 3 - 1.8I+8X10 3 I o (qr))t r 

»j( 1 . 985 xl 0 3 1.366x10* i o (qr))t ( j ) ] 

H = A(<»,t) ^(cjr )t z 
P = A( 4>,t ) [4j.. 57xl0~ 6 I-^qr)) 
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CUTOFF FIELDS: 

H-Wave 
m = 1 
u = 5.77 

p = = j q 

A(<>,t) = A-jO' 55 * 77UJnt 

Complex Fields: 

E = Skk A(*,t) [(1 + (|) 2 )? r + 3d - (f) 2 )^] 

H - AH>,t)(|)t z 

7 = 52.8 A(+,t) [-3(1 + 3-l6(|) 2 )T r + (1 - 3.16(|) 2 )t^l 
p = 3.85 x 10' 7 A(* 1 t)(|) 
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Discrete Frequency: 
m = 0 



u = 1 

8 n = 2.1(45 



"1 

?2 



= 3-415 e j 



ir/4 



= hlUl e J" A 



A( z, t ) = A^e 



j< ui n t-2.445 f) 



Complex Bessel Function 

(pr ) = U^pr) + JV i (pr) 

Conplex Fields: 

E = A(z,t) fj(-3»55 (P^ r ) + 1-068 V 1 (p 1 r))t r 
-(.957 U-^r) + 1.588 V 1 (p 1 r)) 

-**(1.14.90 U 0 ( P;L r) - .666 V o ( P;L r))t z ] 

H = A(z,t) [(6.22xl0~ 3 U-^p-jT) + 10.29xl0“ 3 V^p^))* 
-(17.77xlO“ 3 TJ-^p-jr) + 6.83xlO” 3 V^p.^))*^ 

-j (16.27 x 10” 3 U 0 ( Pl r) + 3.99*10~ 3 V 0 ( P;L r))? z ] 

J = A(z,t) [(.679 (p^r ) + .393 V ] ,(p 1 r))? r 

+j (-1.231 U-^Pjr) + .1724 V 1 (p 1 r})^ 

+j (-1.266 U Q (p 1 r) + .566 V 0 (p 1 r))‘? z ] 

P = A(z,t) [jxlO' 9 (-1.702 U Q ( Pl r) + 2.26 V 0 ( P;L r)] 
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Discrete Frequency: 
tn = 0 

u = 4* 773 
Q n = 1 * 663 

P 1 = 2.4o5/a 
P 2 = 8 . 654/a 

A( 2 ,t, = VJ^ 773 '^- 1 ' 663 S> 

Conplex Fields: 

E = A(z,t) [j(.2l|6 J x (p x r ) = 3.15 J 1 (p 2 r))t r 
+(".337 J 1 (P 1 r) + .843 ^(pgrjjt^ 
+(J 0 ( p l*) + *866 J 0 (P 2 r ))^ z 3 
H = 10~ 3 xA(z,t) [(.312 J x (p ir ) - .594 J^p-jT))^ 

+j (1.568 J 1 (p-jr ) + .298 ^(pgr))^ 

+j (-.451 J^P-jr) + 3.19 J 0 (P 2 r))t z ] 

J = A(z,t) [(.0096 J^p-jT) - .807 ^(pgr))^ 

+j (.0533 J^p^r) + .393 ^(Pgr))^ 

=j (.1779 J Q (p-lX- ) + .1183 J 0 (P 2 r))t z ] 
p = -jlO” 9 A(z,t) [.1883 J^P-jr) + 5.01 J 0 (p 2 r)] 
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